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MnAu 2 is a spin-spiral material with in-plane ferromagnetic Mn layers that form a screw-type pat¬ 
tern around a tetragonal c axis. The spiral angle 9 was shown using neutron diffraction experiments 
to decrease with pressure, and in later studies it was found to suffer a collapse to a ferromagnetic 
state above a critical pressure, although the two separate experiments did not agree on whether this 
phase transition is first or second order. To resolve this contradiction, we use density functional 
theory calculations to investigate the spiral state as a function of pressure, charge doping, and also 
electronic correlations via a Hubbard-like U. We fit the results to the one-dimensional J 1 — J 2 — J 3 — J 4 
Heisenberg model, which predicts either a first- or second-order spiral-to-ferromagnetic phase tran¬ 
sition for different regions of parameter space. At ambient pressure, MnAu 2 sits close in parameter 
space to a dividing line separating first- and second-order transitions, and a combination of pressure 
and electron doping shifts the system from the first-order region into the second-order region. Our 
findings demonstrate that the contradiction in pressure experiments regarding the kind of phase 
transition are likely due to variations in sample quality. Our results also suggest that MnAu 2 is 
amenable to engineering via chemical doping and to controlling 6 using pressure and gate voltages, 
which holds potential for integration in spintronic devices. 


I. INTRODUCTION 

The tetragonal material MnAu 2 is one of the oldest 
known spin-spiral materials m, with a Neel temper¬ 
ature of Tn = 363 K [T] and a local Mn moment of 
Em]. The magnetic ground state consists of in¬ 
plane ferromagnetic Mn layers with a screw-type pattern 
around the tetragonal c axis, which can be described us¬ 
ing the spiral angle 9 [this is equivalent to the wave vec¬ 
tor q = (0,0,52)5 where 9 — qzc/2]. Neutron diffraction 
measurements EllSIS] find that 9 has a weak dependence 
on temperature, increasing slightly with increasing tem¬ 
perature [0(77 K) = 47° and 0(295 K) = 51°[. The spiral 
state is known to collapse to a metamagnetic fan-like con¬ 
figuration when placed in a ^ 10 kOe magnetic field at 
room temperature [1112], which led to a revival of inter¬ 
est in MnAu 2 when it was found that this gives rise to a 
giant magnetoresistive effect j^. 

There are several microscopic mechanisms that can in¬ 
duce spin spirals. In itinerant systems, both Fermi sur¬ 
face nesting |7] and the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interaction |5HTn] can lead to spiral formation. 
The Dzyaloshinsky-Moriya (DM) interaction |lll [12], 
which occurs in materials without an inversion center 
(MnAu 2 has an inversion center, which rules it out as a 
possible mechanism), can also induce spirals and moment 
canting, but due to the interaction’s relativistic origin it 
is more important in materials with heavy elements. It 
can still be competitive in lighter transition metals, such 
as MnSi [THUTH] . but the spirals will have a long wave¬ 
length due to weak relativistic effects. Finally, frustra¬ 
tion, which can be geometric or magnetic, can encourage 
noncollinearity and lead to the formation of spirals. 

The spin spirals in MnAu 2 are understood to be a 
consequence of magnetic frustration, which is commonly 
modeled using a one-dimensional Heisenberg model with 
first and second neighbor interplanar couplings (the Ji — 
J 2 model) |161IT7| . In this model, the spiral state is stable 


when Ji yf 0, J 2 > 0, and jJij < 4 |J 2 | [11]. The valid¬ 
ity of applying this model to MnAu 2 was confirmed us¬ 
ing density functional theory (DFT) calculations [T8l[T9|. 
The origin of the magnetic frustration was traced to a 
competition between two exchange mechanisms, nearest- 
neighbor superexchange and a transferred RKKY-like in¬ 
teraction, and electronic correlations were found to play 
an essential role in suppressing the RKKY-like interac¬ 
tion, which is necessary to satisfy the jJij < 4 |J 2 | in¬ 
equality m- 

The spiral angle 0 is sensitive to pressure, with neu¬ 
tron diffraction measurements determining that 9{P) de¬ 
creases with applied pressure [ 0(0) = 50.7° (47.0°) de¬ 
creasing to 0(8.83 kbar) = 41.8° (40.5°) at a tempera¬ 
ture of 295 K (77 K)[ |S|. Measurements of the critical 
magnetic field He for the spiral-to-ferromagnetic transi¬ 
tion as a function of pressure, when extrapolated to zero 
external field, implied that the spiral state should col¬ 
lapse to ferromagnetism at ~ 12 kbar [20]. The pressure- 
induced spiral-to-ferromagnetic transition was confirmed 
by subsequent experiments using inductance measure¬ 
ments ini, electrical resistivity measurements, and mea¬ 
surements of He |H] , although the two reports disagreed 
on the order of the phase transition. In Ref. HTJ the 
transition is of the second kind and occurs over a pres¬ 
sure range of 12-20 kbar, while in Ref. |21]it is of the first 
kind and occurs at 13 kbar. This contradiction was never 
resolved and remained an open question. 

The Ji — J 2 Heisenberg model, which can explain the 
stability of the spiral state, predicts that the transition 
must be of the second kind. However, general investiga¬ 
tions of the classical one-dimensional Heisenberg model 
show that including the third nearest-neighbor term J 3 
leads to regions of parameter space with first order phase 
transitions |23[[^ . Further neighbors, such as J 4 , change 
the area of parameter space where first-order transitions 
occur 125 . It follows then that if a system were to sit near 
the boundary in parameter space separating first- and 
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second-order phase transitions, then variations in mate¬ 
rial quality could shift the system towards one region or 
the other. 

We argue that this is the case for MnAu 2 . We support 
this claim using DFT calculations to simulate the effects 
of pressure, charge doping, and variation in the strength 
of electronic correlations via a Hubbard-like U. We cal¬ 
culate the magnetic energies and fit them with a clas¬ 
sical, one-dimensional Heisenberg model with coupling 
between the first four nearest-neighbor Mn planes. We 
confirm that ideal, stoichiometric MnAu 2 sits close to a 
boundary separating first- and second-order phase tran¬ 
sitions, and that varying the pressure and doping shifts 
MnAu 2 across this critical boundary. This suggests that 
the two pressure experiments are not in conflict, but in¬ 
stead the MnAu 2 samples from each are located in two 
separate portions of phase space. The spiral angle 9 is 
also found to be responsive to applied pressure and dop¬ 
ing, which holds promise for integration of MnAu 2 into 
spintronic devices. 


II. COMPUTATIONAL DETAILS 

We employed DFT to solve the electronic structure 
of MnAu 2 and calculate the total energy of magnetic 
configurations. To perform these calculations, we used 
the all-electron code elk |^, which is an implemen¬ 
tation of the full-potential linear augmented planewave 
method, and also projector augmented wave potentials 
implemented in the code VASP |26l|27]. We used the lo¬ 
cal spin-density approximation (LSDA) for our cal¬ 
culations, and for correlation effects we used the DFT-|-t7 
method in the fully localized limit |29| . in which we in¬ 
troduce a Hubbard-like U on the 3d orbitals of the Mn 
atoms. We used two values ofU, U = 3.12 eV and 4.7 eV, 
and J was set to 0.7 eV. 

MnAu 2 has a tetragonal crystal structure with 
lA/mmm space group symmetry. The experimental lat¬ 
tice parameters are a = 3.370 A and c/a = 2.599, yielding 
a volume of 49.741 A^/formula unit, and the Wyckoff po¬ 
sitions for Mn and Au are 2a and 4e, respectively, with 
the internal experimental height for Au zau = 0.34 (frac¬ 
tional coordinates) |301 I3T] . To simulate pressure, we 
began with the MnAu 2 unit cell from experiment with 
tetragonal symmetry and varied the lattice parameters 
a and c to set the volume |32]. We then fixed the vol¬ 
ume and performed structural relaxations in VASP to op¬ 
timize the c/a ratio and the internal parameter zau- Af¬ 
ter the relaxation, we performed total energy calculations 
in VASP with the tetrahedron method and calculated the 
pressure using the formula P = —dE/dV. We then im¬ 
ported the relaxed structures into elk for calculating the 
magnetic energies. 

We employed two methods to simulate charge doping, 
( 1 ) adding electrons directly to the system along with 
a uniform positive charge background (hereafter referred 
to as “charge dosing”), and ( 2 ) using the virtual crystal 


approximation (VGA). The VGA involves replacing the 
Au ions with fictitious ions of fractional charge in order 
to add electrons or holes to the system. We used the 
experimental lattice parameters of MnAu 2 in our dop¬ 
ing calculations, while for calculations with simultaneous 
charge doping and applied pressure we used the relaxed 
cells obtained with VASP. 

We used the spin spiral method implemented in ELK to 
calculate the total energy as a function of the spiral angle. 
This method defines a wave vector q = ( 0 , 0 , 92 ) that is 
applied to a primitive cell to simulate spirals. In our 
previous study of MnAu 2 m, we found good agreement 
between the calculated energies obtained via this method 
and noncollinear calculations of spirals in supercells. In 
this study, our procedure was to calculate the energy as 
a function of the wave vector for the range 0 < 9 ^ < 27r/c 
as we varied the pressure, charge doping, and Hubbard 
U. We then used 6 = qzc/2 to convert the wave vector 
to an angle. 

Gonvergence for the spin spiral calculations was care¬ 
fully checked in the primitive cells, and the following pa¬ 
rameters were necessary to achieve convergence: k-mesh 
= 16 X 16 X 16, nempty = 8 , rgkmax = 8.0, gmaxvr = 
12 , Imaxapw = 8 , and the smearing function width was 
set to 0.0001 Ha. 


III. SUMMARY OF ID HEISENBERG MODEL 

The magnetic interactions in MnAu 2 are typically 
modeled using the Heisenberg model: 

= -nij (1) 

Note that = J^jmp, where Jy is the magnetic cou¬ 
pling and |m| is the moment amplitude. DFT calcu¬ 
lations confirm that Mn is in the high spin state with 
S = 5/2, which makes the classical Heisenberg model 
suitable for studying the magnetic interaction. The mag¬ 
netic moment is also assumed to be constant, which is 
reasonable as the Mn moment amplitudes do not vary by 
much as a function of 9. 

The spiral state of MnAu 2 is completely described by 
the parameter 9. As a result, we can simplify the Heisen¬ 
berg model to one dimension: 

N 

H = const. + ^ J„ cos(n0) (2) 

n—1 

The interplanar coupling determines the spiral state 
stability, and below we discuss the phase diagram for 
the models which include two (Ji — J 2 model), three 
(Ji — J 2 — J 3 model), and four (Ji — J 2 — J 3 — J 4 model) 
nearest-neighbors in Eq. (i). 

The Ji — J 2 model is the simplest possible model with 
a stable spiral state, which occurs when Ji 7 ^ 0 , J 2 > 0 , 
and |Ji| < 4 |J 2 |- The angle 9 is plotted as a function 
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FIG. 1. (Color online) Phase diagrams for the one-dimensional Heisenberg model of Eq. Q. (a) The phase diagram of the 
Ji — J 2 model, (b) The phase diagram of the Ji — J 2 — J 3 model, where the contour plot corresponds to 6. The thick black 
line indicates a first-order phase transition and the thick white line indicates a second-order transition, (c) The phase diagram 
of the Ji — J 2 — J 3 — J 4 model for J 4 /IJ 1 I = —0.0346, plotted in the same manner as panel (b). The white star indicates 
where ideal stoichiometric MnAu 2 with U = 3.12 eV sits in the phase diagram, (d) The critical values of J 2 /I Ji| (contour plot), 
J 3 /IJ 1 I, and J 4 /IJ 1 I for the spiral-to-ferromagnetic phase transition. The red line divides parameter space into two regions, 
where first-order phase transitions occur north of the line, and second-order transitions occur south of it. The red and white 
stars represent where ideal, stoichiometric MnAu 2 sits for U = 3.12 eV and U = 4.7 eV, respectively. 


of J 2 /\Ji\ in Fig. [^a), and the spiral-to-ferromagnetic 
phase transition is of second order. What has not been 
appreciated in previous studies of MnAu 2 is that this 
model cannot explain how the turn angle could collapse 
under pressure in a first-order transition, as measured in 

Ref.m 

Including the third nearest-neighbor coupling yields 
the Ji — J 2 — J 3 model, which can be solved analytically 
(see Supplemental Material [S^) H31|M]. The phase di¬ 
agram of 9 for this model as a function of J 2 /\Ji\ and 


J 3 /IJ 1 I is depicted in Fig. Sb). There are two regions 
in the parameter space, the ferromagnetic region in the 
lower left and the spiral region in the rest of the diagram, 
with the contour plot corresponding to the magnitude of 
9. The phase boundaries separating the two regions are 

m-- 

l>{\Ji\-^J2) 0< J2/|Ji|<2/5 

kAt) J2/\Ji\>2/5 

Importantly, this model contains both first- and second- 
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FIG. 2. (Color online) MnAu 2 lattice parameters as a function 
of pressure. 


order phase transitions, which are distinguished by the 
kind of degeneracy that occurs on the borders in Eq. ^ . 
A first-order phase transition occurs when there is a de¬ 
generacy between the ferromagnetic solution (9 — 0) and 
a finite angle (0 > 0 ), while the second order phase transi¬ 
tion corresponds to a smooth, continuous connection be¬ 
tween ferromagnetic and spin-spiral regions, and is found 
by taking a second-order Maclaurin series expansion of 
Eq. (§ and equating it with the ferromagnetic solution 
in the limit 0 —> 0. Second-order transitions correspond 
with the 0 < J 2 /\Ji\ < 2/5 result and are depicted as 
the thick white line in Fig.Wb), while first order transi¬ 
tions correspond to the J 2 /\Ji \ >2/5 result and are de¬ 
picted as the thick black line. From this, it is clear that 
a small ferromagnetic coupling between third-neighbor 
planes (relative to | Ji |) can change the order of the phase 
transition. 

The analysis is more complicated for four or more 
neighbor couplings since an analytic solution for 9 is 
not available. Despite this, some properties can still 
be analytically solved (see Supplemental Material [33]) 
|24| . The phase boundary for second order transitions 
is J 2 = (I Ji| — 9 J 3 — I 6 J 4 ) /4. The criterion for a first 
order phase transition is 




l-k64J4/|Ji| 
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(4) 


Note that the expression for the second order phase 
boundary and Eq. ([^ are valid when J3, J4 > 0. If one 
of the parameters is negative, then the expressions still 
hold as long as the antiferromagnetic couplings domi¬ 
nate HH, such as J 3 > 0 , J 4 < 0 , and IJ 3 I IJ 4 I. 
What’s clear is that a relatively small fourth neighbor 
coupling can have a dramatic effect on the location of 
the first order and second order phase boundaries. To il¬ 
lustrate this, we numerically calculate the phase diagram 
for the Ji — J2 — J3 — Ji model with J 4 /I Ji| = —0.0346, 


which we plot in Fig.j^c). Comparing panels (b) and (c), 
we clearly see that the length of the second order phase 
boundary (white line) has changed, despite the relative 
weakness of the fourth neighbor coupling. We also cal¬ 
culate the exchange parameters for ideal, stoichiometric 
MnAu 2 with U = 3.12 eV, which are J 2 /\Ji\ = 0.3840, 
>^3/I>4i| = —0.0043, Ji/\Ji\ = —0.0346, and indicate 
where MnAu 2 sits with the white circle. The material 
sits very close to the first order transition, and a small 
variation in J 3 /IJ 1 I would drive the system towards fer¬ 
romagnetism. 

In Fig. m we plot a contour map of the critical 
values oi J2/\Ji\, J 3 /IJ 1 I, and J 4 /IJ 1 I for the spiral-to- 
ferromagnetic transition in the one-dimensional Ji — J 2 — 
J 3 — J 4 model. We only consider J 2 > 0 in this figure. 
The contour plot shows the surface in (J 2 /IJ 1 I, J 3 /IJ 1 I, 
J 4 /1 Ji I) space that separates ferromagnetic and spiral or¬ 
der. For example, if you calculate the exchange param¬ 
eters for MnAu 2 and those parameters place the system 
“inside” the surface (J 2 < Jf), then the ground state 
magnetic order is ferromagnetic. If instead the param¬ 
eters place MnAu 2 on or outside the surface (J 2 > 
then the magnetic order will be a spiral state. Note that 
the light yellow region corresponding to = 0 means 
that the system has spiral order for all J 2 > 0. Finally, if 
you can tune the parameters, such as by applying pres¬ 
sure or adding holes/electrons, then the system can cross 
the surface and undergo a phase transition. 

As we saw in panels (b) and (c) of Fig. phase tran¬ 
sitions can either be first or second order. The param¬ 
eters (J 3 /I Ji|, J 4 /I Ji|) control the kind of phase transi¬ 
tion that occurs when you cross the surface depicted in 
Fig.Qd). The red line separates the first- and second- 
order transition regions. If you cross the surface while 
north of the line, then the phase transition is second 
order; if you cross the surface while south of the line, 
then the phase transition is first order. Having estab¬ 
lished this, it is now instructive to indicate where ideal, 
stoichiometric MnAu 2 sits in (J 3 /IJ 1 I, J 4 /I Ji|) space for 
U = 3.12 eV (red star) and U = 4.7 eV (white star). For 
both systems, J 2 > J|, so they have spiral order. What is 
striking here is how close MnAu 2 sits to the dividing line, 
especially when U = 3.12 eV. This shows that changes in 
the exchange parameters induced through pressure, dop¬ 
ing, or impurities could place the system north of the red 
line. Since the system sits close to the dividing line, this 
could explain the discrepancy in experimental pressure 
studies discussed in Sec. H) 

One final point to note is that this analysis assumed 
that the local moments always have an in-plane orien¬ 
tation. Neutron diffraction experiments confirm that 
spin spirals in MnAu 2 are oriented in-plane | 2 |, and our 
calculations of the magnetocrystalline anisotropy energy 
for MnAu 2 with U = 3.12 eV give if(OOl) — i?(100) = 
0.35 meV |51|, indicating the in-plane direction is the 
preferred axis. Yet, it is less clear what would happen 
on the first order transition boundary, where zero and 
finite angle configurations are degenerate. We checked 
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P (kbar) dosing (elec/Au) VGA (elec/Au) 


FIG. 3. (Color online) The dependence of the exchange parameters J 2 /I Ji|, J 3 /I Ji|, and J 4 /I Ji|, see legend, as a function of 
pressure, doping, and U. In the charge dosing and VGA panels, positive values on the horizontal axis correspond to electron 
doping and negative values to hole doping, (a) Pressure dependence with [/ = 4.7 eV. (b) Charge dosing dependence with 
U = 4.7 eV. (c) VGA dependence with U = 4.7 eV. (d) Dependence on Hubbard U. (e) Pressure dependence with U = 3.12 eV. 
(f) Charge dosing dependence with U = 3.12 eV. (g) VGA dependence with U = 3.12 eV. 


this possibility by considering conical solutions. We 
added a canted angle ip and a uniaxial magnetocrystalline 
anisotropy term to Eq. (|^, 


E{6, p) = const + ZIsin^((p) + 


^ 1 

n—1 


2 sin^((/?) sin^ 



(5) 


In this expression, the angles 9 and p remain decoupled 
and the solutions for 9 are unaffected by canting and sim¬ 
ple uniaxial magnetocrystalline anisotropy. If we now re¬ 
strict ourselves to the first-order phase transition bound¬ 
ary, the energy difference between the conical spirals and 
the ferromagnetic configuration is 

E{9, p) - £1(0, p') = D [sin2((/?) - sin^((p')] ( 6 ) 

Therefore in the absence of magnetocrystalline 
anisotropy there is a continuous degeneracy of the 
canting angle. This degeneracy is lifted with the 
magnetocrystalline anisotropy term, although we should 
note that the magnetocrystalline anisotropy parameter 
is small, D — 0.35 meV ^ 4 K, and since the pressure 
experiments 01111122 ] were performed at temperatures 
of 77 K and at room temperature, it is possible for the 
ferromagnetic orientation to align along the c axis. It 
could also be used as a continuous pathway to make 9 
drop to zero (at finite 9 the moments can align parallel 
to c while on the first-order boundary, and then as you 
move away from the boundary they can drop back to 
being in-plane with 0 = 0). This means that the phase 
transition may not be truly first order. However, it 
should be noted that such a transition entails a rapid 
change in the canting angle, which would be difficult to 


distinguish from a first-order transition in experiment. 
For the rest of the discussion, we will continue to refer 
to this as a first-order transition with the caveat that it 
may correspond to a rapid change in the canting angle. 
Regardless, it remains distinct from the second-order 
transition, where 9 continuously decreases to zero. 


IV. DFT RESULTS AND DISCUSSION 


Having established in Sec. HI that MnAu 2 sits close to 
a dividing line that separates first and second order phase 
transitions, we now use DFT to perform a systematic set 
of spin-spiral calculations to simulate the effects of pres¬ 
sure and doping, and the effect of varying U. This will 
reveal how applying pressure and doping makes MnAu 2 
evolve in parameter space. Each set of spin spiral calcu¬ 
lations results in a E{9) curve, and all curves are similar 
to what was obtained in our previous report |19| . Each 
curve was fitted using the Ji — J 2 — J 3 — Ja model, and 
from these we extract the exchange constants (see Sup¬ 
plemental Material m)- 

It should be noted that these calculations and their 
subsequent analysis are all done at zero temperature, a 
choice that we should justify. Measurements of the sat¬ 
uration magnetization at low temperatures (T < 10 K) 
and at room temperature |35l 136] , along with measure¬ 
ments of the susceptibility and the critical magnetic field 
He |35| . show a very weak dependence on temperature. 
As discussed in Sec. |T] the temperature dependence of 
the spiral angle as a function of pressure, obtained using 
neutron diffraction measurements at 77 K and at room 
temperature, is also weak 00. The weak temperature 
dependence of these material parameters suggests that 
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Pressure Ji J 2 J 3 J 4 Doping Type Ji J 2 J 3 J 4 

(kbar) (meV) (elec/Au) (meV) 

U = 3.12 eV 


0.000 

-20.1394 

9.8339 -0.9823 -0.7333 

-0.050 

12.977 

-22.3635 

9.9571 -0.7869 -0.7692 


38.109 

-26.5047 

10.1130 -0.4949 -0.7711 

-0.025 

67.892 

-31.1437 

10.2032 -0.2159 -0.6435 


88.298 

-34.1544 

10.3499 -0.1301 -0.6544 

0.000 

108.353 

-37.0067 

10.4056 -0.0440 -0.5633 


169.709 

-44.9705 

10.4817 -0.3415 -0.5701 

0.025 




0.050 


Dosing -12.5109 6.7832 -0.5134 -0.6678 
VGA -14.1347 7.1857 -0.5786 -0.6446 
Dosing -16.2810 7.1840 -0.2603 -0.6489 
VGA -17.1181 7.3614 -0.3046 -0.6662 
Dosing -20.2505 7.7762 -0.0861 -0.7000 
VGA -20.2505 7.7762 -0.0861 -0.7000 
Dosing -24.5855 8.2909 0.3200 -0.7065 
VGA -23.5647 8.1775 0.3180 -0.6913 
Dosing -29.1444 8.4463 0.9843 -0.5759 
VGA -26.9893 8.2059 0.8370 -0.6678 


U = 4.7 eV 


0.000 -5.5435 

67.892 -13.4226 
169.709 -24.2692 
200.000 -27.4309 
240.000 -31.4087 
317.691 -39.3228 


6.4037 -0.6415 
7.7049 -0.2128 
9.1221 0.0078 
9.6846 0.0290 
10.4061 0.1081 
11.3247 0.2179 


-0.9926 

-0.050 

-1.0301 


-1.0323 

-0.025 

-1.0575 


-1.1649 

-0.010 

-1.1539 

-0.005 


0.000 


0.005 


0.010 


0.025 


0.050 


Dosing 

VGA 

Dosing 

VGA 

Dosing 

Dosing 

Dosing 

VGA 

Dosing 

Dosing 

Dosing 

VGA 

Dosing 

VGA 


-0.6574 

-1.6202 

-3.1144 

-3.6679 

-4.8464 

-5.4763 

-6.1229 

-6.1229 

-6.7465 

-7.4455 

-9.6007 

-8.8932 

-13.4470 

-11.7329 


3.6749 

4.0940 

4.3517 

4.5317 

4.7346 

4.9005 

5.0350 

5.0350 

5.2203 

5.4008 

5.9304 

5.7378 

6.5046 

6.2087 


-0.6583 

-0.7241 

-0.9235 

-0.9799 

-0.9704 

-0.9132 

-0.9103 

-0.9103 

-0.8760 

-0.8414 

-0.7010 

-0.6957 

-0.4418 

-0.4044 


-0.7146 

-0.7442 

-0.3170 

-0.3207 

-0.1601 

-0.0754 

-0.0480 

-0.0480 

0.0486 

0.1156 

0.3580 

0.3137 

0.6154 

0.6055 


U 

(eV) 

Ji 

J 2 J 3 Ji 

(meV) 

1.7 

-37.7314 

9.0945 -0.2716 -0.2806 

2.2 

-31.4231 

8.9775 0.0325 -0.4564 

2.7 

-25.0395 

8.5043 0.0267 -0.6542 

3.2 

-19.4161 

7.6273 -0.0997 -0.7022 

3.7 

-14.3696 

6.7687 -0.0028 -0.7722 

4.7 

-6.1229 

5.0350 -0.0480 -0.9103 


TABLE I. Values of the exchange parameters as a function of pressure, charge dosing, the VGA, and U. Pressure and doping 
results for U = 3.12 eV and U = 4.7 eV are reported in the first two sections of the table, where the first five columns are the 
pressnre results and the last six columns are the doping results. The “type” column indicates whether the doping was performed 
using charge dosing or the VGA. The parameters as a function of U are reported in the last section of the table. 


the internal energy determines the important qualitative 
features of the magnetic phase diagram of MnAu 2 , and 
that including thermal fluctuations will result in only mi¬ 
nor changes to the area of the spiral and ferromagnetic 
regions and the length of the first and second order phase 
boundaries. Therefore, it is reasonable to expect that our 
conclusions about the spiral-to-ferromagnetic phase tran¬ 
sition obtained at zero temperature hold at finite temper¬ 
atures. 

For the pressure calculations, we constrained the vol¬ 
ume of our cell and relaxed the structure as described in 
Sec. 1^ and then calculated the pressure using the for¬ 
mula P = —dE/dV (see Supplemental Material |S1]). 
Using the LSDA functional during relaxation results in 
over binding, and so here the equilibrium volume is pre¬ 
dicted to be 6% smaller than experiment. To account for 


this, we shift our pressures such that P = 0 corresponds 
to the experimental volume, and the lattice and internal 
parameters as a function of pressure are shown in Fig. 

As the pressure is increased, the parameters decrease in 
a smooth way, which contrasts with the experimental ob¬ 
servation in Ref. which seemed to indicate a sharp drop 
in the c parameter for P > 6.5 kbar m- 

The fitted exchange parameters as functions of pres¬ 
sure, charge dosing, the VGA, and U are reported in Ta¬ 
ble m We find that Ji becomes more ferromagnetic with 
increasing pressure, electron doping, and decreasing U, 
and that it is also the most sensitive to these parameters. 
At 28 kbar of pressure, the absolute value of Ji grows by 
3.0 meV - 35 K for U = 3.12 eV and 4.5 meV - 52 K 
for U = 4.7 eV. In Ref. 1^ the critical transition temper¬ 
ature (Neel or Curie, depending on the magnetic state 
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VGA (elec/Au) U (eV) 

FIG. 4. (Color online) The spiral angle 0 as a function of (a) pressure, (b) charge dosing, (c) the VGA, and (d) U. Panels (a), 
(b), and (c) include results for U = 3.12 eV and 4.7 eV, see legend. The angle 6 is obtained by minimizing the Ji — J 2 — Ja — J 4 
model. 


at the given pressure) increases by ~ 40 K going from 
zero pressure to 28 kbar, so these changes in Ji are of 
the correct order. 

To evaluate the spiral-to-ferromagnetic phase transi¬ 
tion in the context of the phase diagrams in Fig. we 
plot the scaled exchange parameters, ji = Ji/\Ji\, in 
Fig. [31 We find that in all cases j 2 is always antifer¬ 
romagnetic and is the most sensitive to pressure, doping, 
and variation in U. J 2 decreases with increasing pressure, 
electron doping, and decreasing U. The exchange cou¬ 
plings js and are closer to zero, and do not respond to 
pressure, doping, and U in the same way as j 2 - Applying 
pressure suppresses ja, and for {7 = 4.7 eV the coupling 
becomes slightly antiferromagnetic for P > 170 kbar. 
Charge dosing and VGA doping have a stronger effect, 
where J3 becomes more ferromagnetic with hole doping 
(for U = 4.7 eV, the coupling becomes enhanced at the 
large hole doping of -0.05 electrons per Au) and turns an¬ 
tiferromagnetic with electron doping. However, is in¬ 
sensitive to U, staying close to zero for all considered val¬ 
ues. j4, unlike J 3 , remains ferromagnetic in all cases, and 
often |j 4 | > Ijaj. Pressure affects j 4 in a similar way as js, 
and for U = 4.7 eV so does doping, but for U = 3.12 eV, 
j 4 is less sensitive to doping than J 3 . Finally, decreasing 
U suppresses J4, going from moderate ferromagnetic cou¬ 
pling at U = 4.7 eV to being nearly zero at U = 1.2 eV. 
This indicates that measuring the fourth neighbor inter- 
planar coupling in MnAu 2 would yield information about 
the strength of the electronic correlations in this mate¬ 
rial. 

Next we plot 0 as a function of pressure, charge dosing. 


VGA doping, and U in Fig. [^ We obtain 9 by minimiz¬ 
ing Eq. (j^ with the fitted exchange parameters [3S]. We 
find that 9 is quite sensitive to pressure, doping, and U. 
As expected, applying pressure decreases 9 in agreement 
with Ref. Ul and at a critical pressure this induces a phase 
transition, which is in the range 12.977 < P < 38.109 
kbar for U = 3.12 eV and 169.709 < P < 200 kbar for 
U = 4.7 eV. This indicates that measurements of the 
critical pressure provides information about the strength 
of the electronic correlations. In addition, for both charge 
dosing and VGA doping the rate of change of 9 is inde¬ 
pendent of U, and 9 decreases with electron doping and 
increases with hole doping. Finally, in Fig. j^d) we find 
that 9 decreases with decreasing U, and eventually col¬ 
lapses to a ferromagnetic state. 

The evolution of the magnetic coupling of MnAu 2 
when subjected to pressure, doping, and variation in 
U can be further visualized using the contour plot in 
Fig. gd). We take the exchange parameters J 3 and J 4 
and plot them as parametric curves on top of the con¬ 
tour plot, which captures how they vary with pressure, 
doping, and U. For each point, we also compare the 
fitted j 2 with the critical that marks the spiral-to- 
ferromagnetic phase transition. If j 2 < j^, we plot the 
point as a circle, indicating that the system is ferromag¬ 
netic, and if j 2 > we plot the point as a star, which 
represents spiral order. This procedure yields Fig.j^ and 
on each panel we give a reference point, which for panel 
(a) corresponds to the relaxed cell at experimental vol¬ 
ume (P = 0 ), for panels (b) and (c) corresponds to the 
experimental lattice parameters with zero doping, and 



























FIG. 5. (Color online) Parametric plots of the MnAu 2 ex¬ 
change parameters as a function of pressure, doping, and U. 
The parametric curves are visualized on top of the contour 
plot from Fig. m- The plot marker shapes indicate the 
magnetic order of MnAu 2 , with stars corresponding to spiral 
order and circles to ferromagnetic order. See legend for the 
values of U. (a) Parametric curve as a function of pressure. 
The reference point is the P — 0 unit cell, see Fig.j^for lattice 
parameters. The arrow indicates the direction of increasing 
pressure, (b) Parametric curve as a function of charge dosing. 
The reference point is a unit cell with experimental lattice 
parameters and no doping. The arrows indicate the direction 
for increasing electron/hole doping, (c) Parametric curve as 
a function of the VGA. The reference point is a unit cell with 
experimental lattice parameters and no doping. The arrows 
indicate the direction for increasing electron/hole doping, (d) 
Parametric curve as a function of U. The reference point is a 
unit cell with experimental parameters and U = 4.7 eV. The 
arrow indicates the direction of decreasing U. 


The parametric path for MnAu 2 under pressure in 
Fig. E a) shows that the spiral-to-ferromagnetic transi¬ 
tion is first order for both values of U. So, for ideal, 
stoichiometric MnAu 2 samples, we expect the pressure- 
induced phase transition to be of first order, in agree¬ 
ment with Ref. [22] Decreasing U at ambient pressure 
also induces a first order phase transition. In contrast, 
for our considered doping levels MnAu 2 does not collapse 
to the ferromagnetic state. Doping with holes stabilizes 
the spiral state and drives the system deeper into the 
first-order parameter region, while doping with electrons 
moves MnAu 2 closer to ferromagnetism and also towards 
a second-order phase transition. Also, the stronger the 
correlations, the more dramatic the shifts in parameter 
space with doping. For both U = 3.12 eV and 4.7 eV an 
electron doping of 0.05 electrons per Au places the sys¬ 
tem at the dividing line between first and second order, 
and it would be expected that the application of pressure 
for this system would induce a second order transition. 

We check the prediction of inducing a second order 
phase transition by calculating the pressure dependence 
of MnAu 2 with U = 3.12 eV and a charge dosing level 
of 0.07 electrons per Au. The exchange parameters ob¬ 
tained from these calculations are reported in Table [ll] 
and the scaled exchange parameters, pressure depen¬ 
dence of 0, and a parametric plot of the exchange pa¬ 
rameters are plotted in Fig. As shown in Fig. I^a) the 
initial charge dosing biases j 2 to be less antiferromagnetic 
and for to become more antiferromagnetic. For the 
range of pressures considered, j 2 decreases with increas¬ 
ing pressure while ja and j 4 are unaffected. Fig. ib) 
shows the effect this has on the 6 vs pressure curve, 
where the trend of the charge-dosed curve is consistent 
with a second order phase transition, particularly when 
compared with the first order transition of the undoped 
curve included for comparison. The parametric curve in 
Fig.|^c) shows that the charge-dosed pressure-dependent 
path in exchange parameter space is north of the dividing 
line, confirming the prediction that applying pressure to 
electron doped MnAu 2 can induce a second-order phase 
transition. 

We can use the result that applying pressure to 
electron-doped MnAu 2 induces a second order phase 


for panel (d) corresponds to experimental lattice param¬ 
eters with U = 4.7 eV. The arrows with text indicate the 
direction that the parametric curves evolve, so following 
the curve from the reference point in panel (a) shows 
how the parameters change with increasing pressure, in 
panels (b) and (c) one direction is electron doping and 
the other is hole doping, and for panel (d) U decreases 
along the curve. If, while following the curve, the marker 
symbol changes from a star to a circle, this corresponds 
to a phase transition. If this change in symbols happens 
while south of the red line, then the phase transition is 
first order, and if it happens north of the red line, then 
it is a second order. 


Pressure Ji J 2 J 3 J 4 

(kbar) (meV) 

0.000 -33.5978 10.2441 0.8442 -0.6193 

2.000 -33.9957 10.2223 0.8197 -0.6252 

6.000 -34.8238 10.1037 0.8679 -0.6570 

10.000 -35.6126 10.0322 0.8756 -0.6331 

18.000 -37.0886 9.9025 0.9327 -0.6302 

20.000 -37.4883 9.8482 0.9131 -0.6498 

24.000 -38.1828 9.7781 0.9236 -0.6241 


TABLE II. Values of the fitted exchange parameters for the 
simultaneous pressure and charge dosing calculations. Results 
are for a charge dosing value of 0.07 electrons per Au and 
U = 3.12 eV. 
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FIG. 6. (Color online) Plots showing the results of calculations with simultaneous charge dosing and applied pressure. The 
calculations were performed using U = 3.12 eV and a charge dosing of 0.07 electrons per Au. (a) The exchange parameters 
J 2 /I Ji|, J 3 /I Ji|, and J 4 /I Ji|, see legend, as a function of pressure, (b) The spiral angle 0 as a function of pressure with and 
without charge dosing, see the legend, (c) Parametric plot of the exchange parameters as a function of pressure visualized on 
top of the contour plot from Fig. [^d). The plot markers indicate the type of magnetic order, see the legend. The reference 
point is the relaxed P = 0 unit cell. The arrow indicates the direction of increasing pressure. 


transition to resolve the contradiction in the pressure ex¬ 
periments. In Ref. HU where the authors found a sec¬ 
ond order phase transition for MnAu 2 , they reported 
that there were trace amounts of MnAu in their sam¬ 
ple, such that the system was 34.3% Mn and 65.7% Au. 
This excess of Mn in their sample can lead to an effec¬ 
tive electron doping. Let us assume that the excess Mn 
atoms are also in the 2 -|- ionic charge state, and that the 
extra electrons they donate are spread homogeneously 
over the system. This would mean that each unit cell 
of MnAu 2 would have 2 electrons per Mn • (34.3/33.3) = 
2.06 electrons per Mn, or an excess of 0.03 electrons per 
Au. Assuming that U = 3.12 eV 133 and comparing 
with the electron dosing curve in Fig.j^b), an extra 0.03 
electrons per Au would shift the sample near the second 
order transition region of the parameter space, which as 
we’ve seen can lead to a second order phase transition 
upon application of pressure. While the MnAu domains 
are not likely to uniformly dope the material, we see how 
the variations in sample quality can influence the kind of 
phase transition. This both stresses the importance of 
controlling for impurities in MnAu 2 , and that the mag¬ 
netic properties are also amenable to tuning. 


V. CONCLUSIONS 

We studied how the magnetic exchange in MnAu 2 
changes upon application of pressure and doping using 
a series of DFT calculations, as well as varying the elec¬ 
tronic correlations with a Hubbard-like [/. We found 
that the proper model for these interactions is the one¬ 
dimensional Heisenberg model with coupling between the 
first four Mn planes (the Ji — J 2 — J 3 — J 4 model), 
which contrasts with previous studies that only consid¬ 
ered interactions between the first two planes. Within 


this model, transitions between the ferromagnetic and 
spiral phases can be either first or second order, and J 3 
and J 4 control the order of the phase transition. Our 
analysis finds that ideal, stoichiometric MnAu 2 sits near 
a dividing line in (J 3 , J 4 ) space that separates first- and 
second-order phase transitions. Applying pressure col¬ 
lapses the spiral state as expected, and for ideal, sto¬ 
ichiometric MnAu 2 the transition is first order. Hole 
doping the material further stabilizes the spiral and the 
tendency towards a first-order transition. In contrast, 
electron doping the material makes the spiral less stable 
and moves the system in (J 3 , J 4 ) space closer to a second- 
order transition. We confirm with additional DFT calcu¬ 
lations that applying pressure to electron-doped MnAu 2 
will induce a second-order phase transition. Therefore, 
impurities that provide an effective electron doping, such 
as excess Mn, can cause the spiral-to-ferromagnetic phase 
transition to become second order under pressure. These 
results indicate that MnAu 2 is tunable, where the kind 
of phase transition can be set and the spiral angle con¬ 
trolled through both pressure and doping, which can in¬ 
clude gate voltages. The ability to control the spiral an¬ 
gle holds promise for integration of MnAu 2 in spintronics 
devices, such as in a spin valve. Attaching a thin film of 
MnAu 2 to the free magnetic layer of a spin valve and then 
applying a gate voltage could allow for a current-free ma¬ 
nipulation of the relative orientation between the free and 
pinned magnetic layers, which would reduce Joule heat¬ 
ing and could help in developing smaller-scale electronics 
devices. 
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SUPPLEMENTARY MATERIALS 

The supplemental materials for this report are 
provided in the IPython notebook format. The note¬ 
book file and the data used for analysis are provided 
at https;//github.com/jkglasbrenner/ipython_ 


nb_collapse_and_control_of_the_mnau2_spin_ 
spiral_state_through_pressure_and_doping For 
a readable HTML version of the notebook, download 
the following file and open it in your web-browser: 
https;//github.com/jkglasbrenner/ipython_ 
nb_collapse_and_control_of_the_mnau2_spin_ 
spiral_state_through_pressure_and_doping/raw/ 
master/ipython_nb_MnAu2_pressure_supplemental_ 
material.html 
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